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Abstract 

s ! \ We propose a simple construction of the non-Hamiltonian dynamical systems possessing an 

invariant measure. These non-Hamiltonian systems are deformations of the Hamiltonian systems 
associated with trivial deformations of the canonical Poisson brackets and time transformation. 
Sometimes these non-Hamiltonian systems can be identified with known nonholonomic systems 
and, therefore, the results are illustrated through the Chaplygin ball and the Veselova system. 

■ 1 Introduction 

We consider simple non-Hamiltonian integrable systems which can be obtained from Hamiltonian 
systems by a proper change of variables, including time. The main examples are one of the most famous 
solvable problems in nonholonomic mechanics, describing rolling of a balanced ball over a horizontal 
surface without slipping, is referred to as the Chaplygin ball [10 and the nonholomic Veselova system 

£h ' Let us start with a system of Hamiltonian equations of motion 

d 

— Zk = Z k = {H,z k \, fc = l,...,m, 

where {., .} being canonical Poisson brackets. If g{z) is a nowhere vanishing smooth integrable function, 
' the time reparameterization of the form 

. dt -> g{z)dt 



changes the invariant probability measure, despite the fact that the orbits remain invariant and the 



ergodicity is preserved. Physically, this reflects the fact that the transformed system evolves at different 
speeds and hence with different residence times along the orbits. 

After time reparameterization one gets the desired simple non-Hamiltonian system possessing an 
invariant measure. Associated with the new time vector field 

X 

H ■ Z = g(z)Z 

C3 

is a so-called conformally Hamiltonian vector field. If the initial Hamiltonian vector field Z is complete, 
then the conformally Hamiltonian vector field Z is a Hamiltonian with respect to another Poisson 
bracket {., .}'. If another technical condition is satisfied, then there is symplectic diffeomorphism 
which maps each orbit to itself and is an equivariant with respect to the flows of the vector fields Z 
and Z H7J. 

In classical mechanics transformations Z — » Z and dt — > gdt are very useful to investigate the 
trajectory equivalent systems, see examples in [TTl [23j EU EE]. In nonholonomic mechanics we have 
to add only one intermediate step to this standard scheme and to deform canonical Poisson bracket 
before the time reparametrization. This modification plays a key role in reducing non-Hamiltonian 
systems to Hamiltonian ones [5] [7] |51 [TO1 IT! } 113 1 IT5 1 115 ] . 

If we have some Hamiltonian vector fields Z^\ . . . , Z^ m \ we can consider some more complicated 
time reparametrization associated with the following vector field 

Z = 9l (z)Z^ + ---+g m (z)Z^. 

Such vector fields are useful for Hamiltonization of some nonholonomic systems [26]. 
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Below we prove that transformation of the non-Hamiltonian vector fields Z — > (Z±, Z2) — > Z relates 
the Veselova system with a simple Hamiltonian system on the sphere and then with the nonholonomic 
Chaplygin ball. It allows us to prove that nonholonomic Veselova and Chaplygin systems are trajectory 
equivalent systems. 



1.1 Simple non-Hamiltonian systems on cotangent bundles 

Let Q be an n-dimensional smooth manifold. Its cotangent bundle T*Q is naturally endowed with the 
Liouville 1-form and symplectic 2-form = d6, whose associated Poisson bivector will be denoted 
with P. In local symplectic coordinates on T*Q 

z = (q,p) = (q 1 ,...,q n ,Pi,---,Pn) 
they have the following local expressions 

n n d d 
9 =pxdqx + ...p n dq n , £1 = ^ dq { A dpi , and P = — A - — , (1.1) 

The 2n volume form 

/1 = O™ ~ dqdp , 

is invariant with respect to Hamiltonian evolution defined by the canonical Poisson bracket {.,.}, the 
Hamilton function H and the corresponding time t 

j t q k = {H,q k }, j t p k ={H,p k }, k = l,...,n. (1.2) 

Let us change the form of these Hamiltonian equations (jl.2l) using trivial deformation of the canonical 
Poisson bracket. Namely, let us substitute scaling momenta 

Pk^g(q)p k k = l,...,n, (1.3) 

into the Liouville and symplectic forms 

-> g = g(q) (pidqx + . . .pndq^ , fi -> Vt g = d9 g . 
The corresponding Poisson bracket 

{lu1j}g = . {luPj}g = 9{q)8ij , {Pi,Pj} 3 = F id , (1.4) 

where 

Fij = djgpi - digpj , d k = d/dq k , 

is a trivial deformation of canonical Poisson bracket in the Poisson-Lichnerowicz cohomology [15] . The 
corresponding volume form 

(i g =n^g- n dqdp (1.5) 
is invariant with respect to a Hamiltonian flow associated with time variable t g defined by 



± z = {H,z} g . (1.6) 



Here and below we suppose that g(q) is a nowhere vanishing smooth function on Q. We consider these 
equations as an intermediate deformation of the initial Hamiltonian equations of motion (|1.2|) . 



Proposition 1 Transformation of momenta p — > g{q)p followed by a change of time 

dt d = g(q)dt g (1.7) 
preserves initial equations of motion on q k and shifts equations on p k 

dq k dq k dp k dp k 

-JT = ^r> ~JT = 17 + Xkqk ' k = l,...,n, (1.8 

dtd dt did dt 
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where 

q k g dp 3 

These non-Hamiltonian equations of motion have invariant volume form 

V d = g 1 - n &q&p, (1.10) 
with measure depending on a nowhere vanishing smooth function g(q). 
Proof: The proof is a straightforward verification of the first half of equations 

A qk = g -i {H , q k}g = g-l ^{p,, q k}g °* = {H, q k } = ^q k , 

3 = 1 J 



and calculation of the remaining part of equations 



dt d 



— p k = g 1 {H 1 p k }g = g 



n ( dH 0H\ 



d -lV^r i dH d Pk . x 

j = l r J 

The Jacoby last multiplier g 1- ™ (|1.10|) of equations of motion (|1.8[) satisfies to the standard equation 



EO i_ n azj _ \ -> a _„ azj _ 
di^ ~ ^dZ 9 dH ~ ' 



because g -n is an invariant measure (jl.5p of the Hamiltonian equations of motion (11.61) . □ 

Functions Xi(z) in (|1.8p may be considered as undetermined Lagrange multipliers of some non- 
holonomic system [7] . On the other hand flow (|1.8|) may be associated with the magnetic geodesic flow 
on the Riemannian manifold [5]. 

Example 1 Equations of motion for nonholonomic LR-systems [29] are deformations of the standard 
Euler equations on the Lie algebra g 

±N i = [u,N i ], j t M=[u,M] + \ i N i , (l.ll) 

in which we shift the second part of initial Euler equations on A,;7V l similar to (jl.8[) . Here M £ q* , 
[., .] means a coadjoint action and functions Xi(M, N l ) are undetermined Lagrange multipliers. 
If the singular orbit of q is symplectomorphic to some cotangent bundle T*Q we can try to identify 
equations (fL"TT|) with (TQ|) . 

Natural generalization of the transformation (11.31) looks like 

n 

Pk^^L ki (q)pi k = l, (1-12) 
i=i 

so that 



fc=l 



Here is a trivial deformation of the initial Poisson bracket associated with a (1,1) torsionless 

tensor field L(q), see [3T]. The generic change of time variable t g — > ta (|1.7I) 

— = gi{H 1 ,z} L + g 2 {H 2 , z} L + ■ ■ ■ + g n {H ni z} L , (1.14) 
aid 
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depends on functions gk(z) and integrals of motion Hi, . . . , H n of the initial Hamiltonian system of 
equations (|1.2j) . 

The Turiel deformation (|1.12[) of the initial Poisson bracket yields the following intermediate 
deformation of the invariant volume form 

fi — > fig = det L(q) dq dp , 

whereas transformation of time (|1.14j) leads to the next change of this form [i g — > /id depending on 
integrals of motion Hi and functions gk(z). 

The corresponding example of nonholonomic systems describing a rolling of dynamically asym- 
metric and balanced ball over an absolutely rough fixed sphere may be found in [26] . 

1.2 Momentum map 

If we suppose that some Lie group G acts on M, = T*Q via symplectomorphisms, we can try to transfer 
non-Hamiltonian equations of motion (ll.8[) to the deformed Euler equations on the Lie algebra 

q* using symplectic reduction and momentum map theory [T]. 

For instance, let J 1 *^ -1 be the cotangent bundle of the sphere S n_1 C K." realized as the 
subbundle of S"" 1 x K", 

r*s n_1 = {x, P e x r"; ( P) x) = 0} , 

where (.,.) means the standard scalar product in R™. Following to p] let us consider the standard 
momentum map 

tp : (p,x) G T*^ 1 - 1 -> (L,x) G e*(n) = ao(n) x R" , 

defined by 

Lij — XiPj XjPi . (1.15) 

This standard mapping relates cotangent bundle T*§™ -1 with a singular orbit of the Lie algebra 
e*(n) = so(n) x R™. Remind that e(n) is an algebra of the Euclidian group E(n) of motions of the 
n-dimensional Euclidean space R™. 

It is well known that mapping (|1.15[) is a Poisson map, which transforms canonical Poisson bracket 
(|L~T|) on T*§"" 1 to the Lie-Poisson bracket on e*(n): 

{Lij,Lki} = SkjLii — SuLkj + SjiLki — SkiLji , 

(1.16) 

{L l j,xk} = 5ikx 3 - 5-jkXi , {xi,Xk} = 0. 

It is easy to prove that mapping (|1.15p is a Poisson map for the deformations of the canonical Poisson 
bracket {., .} g (|1.4p or {., .} L (|1 . 13[) on cotangent bundle of the sphere T*§™ -1 , too. 

Proposition 2 Momentum map (1.15\) transforms the Poisson bracket {.,.} g jl-4\ ) on T*§ n_1 into 
the Poisson bracket 

{Mij,M kl } g = g^SkjMn-SiiMkj + SjiMki-SkiMj^ 

+ Mkl }XSi m Xj - S Jm X l )d m g - Mij 2_J^km x i - $i m Xk)d m g , 
m m 

{Mij,xk} g = g{&ikXj - SjkXij , {xi,Xk} g = 0- (1-17) 

on the whole space e*(n). As above here we put d m = d/dx m - 

Proof is a straightforward verification of all the Poisson bracket properties. □ 

On a singular orbit of e* (n) the Poisson brackets (|1.16[> and (| 1 . 1 7[) are related by a trivial change 
of variables 

M i:j = g{x)Lij , = g- x {x)Mij . (1.18) 

We suppose that on generic orbit these brackets are related by more complicated transformation 
discussed in the next paragraph. 
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Example 2 At n = 3 we can identify an element of so(3) with a vector of the angular momentum 
using well known isomorphism of Lie algebras (so(3), [., .]) and (R 3 , x): 

Li — £ijkLj k . 

Here Sij k is a totally skew-symmetric tensor, [a, b] means commutator in so(3) and a x 6 is a vector 
cross product in K 3 . 

In this vector notation the canonical Lie-Poisson bracket f|l . reads as 

{Li , Lj } = SijkLk , {Li ,xj} — SijkXk , {xi,Xj}=0, (1-19) 

whereas the deformed Poisson bracket (| 1 . 1 7|) has the form 

{Mi , Mj } g = e ijk (g(x)M k + x k ^ M m d m g(x)j , 
{M % ,Xj } = e ijk g(x) x k , {xi , Xj } = . 



(1.20) 



These brackets have the following Casimir elements (x, x), (x,L) and (x,M), respectively. 

Momentum map (|1.15l) is a symplectomorphism of T*§ 2 and a singular orbit of e*(3) defined by 

(x,L) = 0. 

On this singular orbit transformation 

L l =g-\x)M l (1.21) 

relates the deformed Poisson brackets (|1.20p with canonical one (j!.19[) . In the next paragraph we will 
consider similar reduction of the Poisson brackets on the generic orbits at (x, L) ^ 0. 

1.3 Simple non-Hamiltonian systems on semi-direct sum of Lie algebras 

Let g be a semisimple Lie algebra, g* its dual space and g = f + p its Cartan decomposition. According 
to [20 , let us consider an affine Poisson subspace of the Lax matrices 

L(X) ~ aX + £ + sA _1 , I G f , s e p* , (1.22) 

passing through the a <E p*. The corresponding Poisson structure on this subspace is isomorphic to 
the Lie-Poisson structure on the semidirect product f* x p*: 

{%4rn}ef*, {%Sfc}€p*, {Si,Sj} = 0. 

The constant a does not affect kinematics, but is quite useful to produce interesting Hamiltonians for 
which equations of motion have the form of the Euler equations 

j t s = {H,s}=[sM, j t l = {H,i}=[ZM + [sA, (1-23) 

where "angular velocity" u> — ip(£) depends on the choice of the Hamiltonian H, whereas variable 
s€p* plays the role of the "Poisson vector" . 

All the details may be found in the book [20], another construction of the integrable systems on 
a semidirect products of a semisimple compact Lie algebra and a commutative algebra may be found 
in [3J[21[T3], see also references within. 

Using semisimple structure we can suppose that a natural generalization of the transformations 
(|1.3I1.18I) looks like 

£ ij ^g(s)l ij + hi j {s,c), (1.24) 

where g{s) is a nowhere vanishing smooth function on commutative part p*, and hij are functions on 
p* depending on the Casimir elements c = (ci, . . . , c m ). Such transformation allows as to change initial 
Poisson brackets {., .} — > {., .} g on a semidirect product of the Lie algebras similar to (|1.4p . 

Now we can add transformation of the time (jl.7p to the transformation of the variables (|1.24[) in 
order to get non-Hamiltonian deformations of the initial Euler equations of motion (|1.23[) : 

-£-s = g- 1 {H,s} g = [s,u], ■^-e = g- 1 {H,£} g = [e,u] + [s,b] + X k s k . (1.25) 

Functions Afc may be considered as undetermined Lagrange multipliers associated with some constraints 
added to the initial equations (|1.23j) . 
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Example 3 Let us continue to consider the Lie algebra e*(3). On its generic orbit the counterpart of 
transformation (ll.3ll.2ip has the following form 

Li = g- 1 (x)^M l -bh l (x)^ , b=(x,M), (1.26) 

where three functions hi{x) satisfy to one algebraic equation 

3 

cg(x) = 1 + Xjh(X) , <^ (x, L) = c (x, M) , ceR, 
»=1 

and one differential equation 

3 



\h 2 - dihi + hi y^^Xjdjhj - Xidjhj) 
i=l j=ii 



= 0, 

ox k 



In this case mapping (|1 .26[) reduces deformed Poisson bracket (|1.20[) to canonical Poisson brackets 
(|1.19p on generic orbits of e*(3). In Section [3] we present particular solutions of these equations 
associated with the Chaplygin nonholonomic ball. 



2 Integrable non-Hamiltonian systems on the sphere 

Transformations (|1.3|) and (|1.7j) change the form of the equations of motion and preserve the main 
invariant geometric and topological properties of the Hamiltonian system (|1.6j) . For instance, it is clear 
that a change of time (jl.7j) or (|1.14[) preserves integrability. In order to get such integrable Hamiltonian 
and the corresponding non-Hamiltonian systems we can use the Jacobi separation of variables method 
[I], method of the classical r-matrix [20], the argument shift method [3j|4l[T4] etc - 

For instance, in the Jacobi method we can substitute some known variables Ui, . . . ,u n (ellip- 
tic, parabolic, conical, prolate spheroidal or other) and conjugated moments p Uk into the separated 
relations 



<$>i(v,i,p Ui , Hi, . . . , H n ) = , i = l,...,n, so that det 



jtO. (2.1) 



Solving the resulting system of algebraic equations with respect to Hi, ... , H n one gets integrable 
Hamiltonian system, which may be interesting in application. After that we can get integrable non- 
Hamiltonian systems (|1.8[) using a proper change of time. 

Let us consider few integrable Hamiltonian systems (|1.2[) and non-Hamiltonian systems (|1.8I ) 
separable in an orthogonal coordinate system on a sphere §™ _1 . Remind, that all orthogonal separable 
coordinate systems can be viewed as an orthogonal sum of certain basic coordinate systems [16) . All 
these basic systems on § n_1 may be obtained from the elliptic coordinate system u — [ui, . . . , u n -i) 
with parameters ei < < ■ ■ ■ < e„ defined through the equation 

~[ X ~ e ^ lli=i(A-e l ) 

Here x t = 1- Like the elliptic coordinates in R n , the elliptic coordinates on S™ _1 are also orthogonal 
and only locally defined. They take values in the intervals 



e\ < Ui < e 2 < u 2 < ■ ■ ■ < u n -i < e„ . 

The coordinates and the parameters can be subjected to a simultaneous linear transformation X{ — > 
axi + b, — > aei + b, so it is always possible to choose ei = and e n = 1. The coordinates can be 
degenerated by letting some, but not all, parameters coincide. 

Let us substitute elliptic coordinates u and conjugated momenta p u into the uniform separated 
relations of the Stackel form 

n 

f( Uj ) H( Uj - ei)p 2 u . = u]- 2 H n _i + ... + u 3 H 2 + Hi + V( Uj ), j = l,...,n-l, (2.3) 
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where rational functions /(it) and V(u) define the Stackel metric and the potential [22] [23] , respectively. 
Solving resulting equations with respect to H±, . . . , if n _i one gets the desired tuple of integrals of 
motion on T*S for the Hamiltonian (|1.2j) and non-Hamiltonian equations of motion (11.81) . 

It is easy to see, that generic separated relations (|2.3[) may be obtained from the separated equa- 
tions with f(u) = 1 

n 

Y[( Uj - ei)p 2 u . = u r ;- 2 H n ^ + ...+ u 3 H 2 + H X + V{ Uj ), j = l,...,n-l. (2.4) 

i=l 

using the transformation 

Vu 2 -> \J ' f{u,)p U] , (2.5) 

which changes the so-called Stackel matrix and the corresponding metric (23] [24]. Below we will 
combine this transformation (|2.5p acting on the form of the separated relations with the similar trans- 
formation (|1.3p changing the symplectic structure and with the time reparametrization (|1.7j) . 



Example 4 Separated relations (12.41) with /(it) = 1 yield many well-known integrable two-dimensional 
Hamiltonian systems on a sphere including the Euler top, Neumann system, Braden and Rosochatius 
systems. For instance, if n=2 and V — 0, the solutions i?i,a of the corresponding separated equations 

3 

¥>3(Mj)i^ =-u :j H 2 +H 1 , ip 3 (u) =4jJ(u-ei), J = 1,2, 

coincide with the Euler integrals of motion 

ifi = (L,AL), H 2 = {L,L), A = diag{e x , e 2 , e 3 ) , 
where L is the angular momentum vector with the entries 

Li = _^V^2 — ___/'( ei _ Ul )p _( ei _ U2 )p V i = 1,2,3, (2.6) 

1*1 — 1*2 V / 

At arbitrary n and U description of the corresponding integrable systems on § n may be found in [2] [30] • 
Example 5 If we combine the transformations (|1.3j) . (|1.7|) and (|2.5|) at 



{ex + !)■■■ (e„ + 1) ' a v v ' y ( ei + 1) • • • (e„ + 1) ' 

one gets the n-dimensional generalization of the nonholonomic Chaplygin ball. Here e(A) is a rational 
function incoming into the definition (|2.2I) of the elliptic coordinates on the sphere. 

For instance, let us put n = 2 and V = 0. In this case the solutions f/i i2 of the corresponding- 
separated relations 



■ ¥>3(%)pL = -«,-(ff2 - dH x ) + dHx , j = l, 2, 



(ei+l)(e 2 + l)(e 3 + l) 
coincide with the integrals of motion for the nonholonomic Chaplygin ball (13.61) 

H 1 = (M,A g M), H 2 = (M,M), 

if ej = dIJ 1 following to [10]. Here Ij are elements of the inertia tensor, d is some parameter and M 
is an angular momentum 

2g(x) XjXk(ej - e k ) 



Mi 



■ ((ei - ui)p Ul - (ei - u 2 )p u ^j , 



1*1 - 1*2 

so that (x, M) = 0. All the details will be presented in Section [3] 

Nonholonomic systems on the sphere with V ^ and the corresponding Lax matrices are consid- 
ered in [35] [55] . Discussion of the n-dimensional generalizations of the Chaplygin ball may be found 
in [11] [M] . Note that the restrictions on parameters ej in [TT] [14] are related with a special choice of 
the nonholonomic constraint. 

Using other functions f(u) (|2.5p , V(u) (|2.4[) and g(u) 111.311.71) we can get a lot of non-Hamiltonian 
integrable systems on the sphere. For instance, if / = u and g = V/i ' ' ' Tni then we get the so- 
called Veselova system and its generalizations, see Section [4] The main problem is how to find non- 
Hamiltonian systems, which are really applicable in mechanics. 
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3 Nonholonomic Chaplygin ball 



The nonholonomic Chaplygin ball [10 is a dynamically balanced 3-dimcnsional ball that rolls on a 
horizontal table without slipping or sliding. 'Dynamically balanced' means that the geometric center 
coincides with the center of mass, but the mass distribution is not assumed to be homogeneous. 
Because of the roughness of the table this ball cannot slip, but it can turn about the vertical axis 
without violating the constraints. 

Let ui and v be the angular velocity and velocity of the center of mass of the ball expressed in the 
so-called body frame. Its mass, inertia tensor and radius will be denoted by m, I and a respectively. 
This reference frame is firmly attached to the body, its origin is located at the center of mass of the 
body, and its axes coincide with the principal inertia axes of the body. The inertia tensor of the body 
in this frame is diagonal I = diag^, Ii, 7 3 ). 

Since slipping at the contact points is absent, its velocity vanishes and we have the following 
nonholonomic constraint 

» + uxr = 0. (3.1) 

Here r is the vector joining the center of mass with the contact point expressed in the body frame and 
x means the vector product in R 3 . 

The angular momentum M of the ball with respect to the contact point with the plane is equal 

to 

M= (I + dE)u) - d(x,u)x, d = ma 2 . (3.2) 

Here x is the unit normal vector to the fixed sphere at the contact point so r = —ax, E is the unit 
matrix and (., .) means the standard scalar product in R 3 . 

Using constraint (|3.1[) we can eliminate the Lagrange multiplier A from the initial equations of 
motion 

mv = mv x oj + A , lib — Iuj x oj + Xr , x — x x uj , 
and obtain the reduced equations of motion in the x, M- variables 

M = M x oj , i = ixw. (3.3) 

All the details may be found in jTOl [7] . 

For further use we rewrite the relation (|3.2p in the following equivalent form 

uj = A 9 M = (a + drf 2 A(x®x) A\M , (3.4) 

where 




A = (I + dE) _1 = a 2 , g = yjl - d{x, Ax) . (3.5) 



There are four integrals of motion 

Hi = (M, A g M) , H 2 = (M, M) , H 3 = (x, x) , H A = (x, M) (3.6) 
and an invariant measure 

^ch — 9~ 1( ixdM , g = — d{x, Ax) . (3-7) 



Thus, the six equations (|3.3[) are integrable in quadratures by the Euler-Jacobi theorem. 

According to O [7] , integrals of motion p.6j) are in involution with respect to the Poisson bracket 

r i ( d(M,Ax) \ r , r , , . 

{Mi,Mj} g = Eijk \ gMk x k j , {M i ,x j \ g = E ijk gx k , {Xi,Xj} g = 0, (3.8) 

where e^k is a totally skew-symmetric tensor. This bracket coincides with the Poisson bracket (|1.20l) 
obtained from (|1.3j) by using a standard momentum map. 
Equations of motion f|3.3[) have the form (|1.25|) 



%=9-HH,z}„ # = (3-9) 



8 



After changing the time variable fj 1 . T|) 

dt -> gdt (3.10) 

these non-Hamiltonian equations of motion are reduced to the Hamiltonian equations with respect to 
the deformed Poisson bracket (I3.8[) . 

We can reduce this trivially deformed Poisson bracket to canonical one. 



Proposition 3 The Poisson map 26) ) 

L a = g- 1 (m 1 - 



bxi ( x\ \ \ cxi 



(x,x) \ v 



L 3 = g-' M; 



bx 3 ( x\ + x\ 



(x,x) 



reduces the deformed Poisson bracket \3. 8\) to a canonical Li-Poisson bracket \1.19\) on e*(3). Here 
the Casimir elements are equal to 

(x, L)=c, (x, M) = b , 

and 

v = x\ + x\ — d(x, x){a\x\ + a2x\) . 

The proof is a straightforward calculation of the Poisson brackets (|1.19j) . □ 

Summing up, using the transformation of coordinates including time we can reduce the initial 
non-Hamiltonian dynamical system Q3.3P to a new Hamiltonian dynamical system on e*(3) with a 
canonical Poisson bracket. In the next paragraph we will study some properties of this Hamiltonian 
system. 

3.1 Hamiltonian system on sphere associated with the Chaplygin ball 

The Poisson map 

M ( L CXl \ ( I X ^ 

1 9 V 1 x\ + x\) (x,x) \ v 

M^g(L^^) + J^(l + X l). 



j i 



(x,x) 



bxs ( x\ + x\ 



M 3 = gL 3 + 1 

is an inverse transformation to p. lip . Applying this map to initial integrals of motion (|3.6p one gets 
integrable Hamiltonian system on e*(3). 

Moreover, at c = one gets a particular Poisson map (|3 . 1 1 1) . which relates generic symplectic leaves 
of the deformed Poisson bracket {., .} g (|3.8|) and singular symplectic leaves of canonical Lie-Poisson 
bracket {., .} (fTT9l) 

(x,L) =0, 

which are symplectomorphic to the cotangent bundle T*S 2 of the sphere. 

In this case initial integrals of motion (|3.6|) are reduced to the following integrals of motion on 
T*§ 2 



n 2/ r A T \, 2b 9(i \ t\ T \ , b 2 {a 1 x 1 +a2X2)(l-d(a 1 x 2 1 +a 2 x 2 2 )) 
Hi = g [L, AgL) + —— y(x, AL) — a 3 x 3 L 3 J H — 



(3.13) 



H 2 = g 2 (L,L)+ 2b9X3L3 + b 2 fi + flM+M) 
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In terms of elliptic coordinates ui,2 (|2.2j) on the sphere and the corresponding momenta p Ul . 
these integrals looks like 

H\ = a P 2 ui + Zbp Ul p U2 + cpl 2 + dp Ul + ep U2 + f , 

(3.14) 

H 2 = Ap 2 ai + 2B PuiPu2 + Cp 2 U2 + Dp Ul + Ep U2 + F . 

Here = dl^ 1 and the coefficients of the polynomials in (|3.14l) are analytical functions on u± t 2- 

According to [T5] , the quasi-Stackel integrals of motion (|3.14[) admit partial separation of variables 
in coordinates s± t 2 defined as the roots of the equation 

(B - bs) 2 - (A - as)(C - cs) = . 

In our case we have 



Sl,2 



d(l + 1(1,2)% si- 



Similar to the Clebsch top [T8] , for the Chaplygin ball these variables evolve on Prym variety of some 
non-hyperelliptic curve, in contrast with the Chaplygin variables of separation [101 [5] . 

4 The Veselova system 

The Veselova system [57J describes motion of the rigid body with a fixed point under constraint 

(u,x)=0. (4.1) 

Here u> — AM is the vector of the angular velocity in the body frame, a; is a unit vector which is fixed 
in a space frame and (., .) denotes the standard scalar product in R 3 . 
The corresponding equations of motion are equal to 

M = M x ui + Xx , x^jxw. (4.2) 

Here M is the angular momentum vector in the body frame and A is the diagonal non degenerate 
matrix 

/ fii 
A = [ d 2 

V d 3 

In order to eliminate the undetermined Lagrangian multiplier A we have to differentiate the constraint 
(|4.1[) and find A using equations of motion (|4. 2[) 



_ (AM x M, Ax) 

A — . 

(Ax, x) 

The resulting equations of motion have four integrals of motion 

Hi = (M, AM) , H 2 = (M,M)~(x,x)- 1 (x,M) 2 , H 3 = (x,x), H 4 = (x,uj), 
and an invariant measure 

fj, ves = g dxdM , g(x) = {x,Ax) 1/2 . (4.3) 

If the constraint (|4.1j) holds then these integrals of motion are in involution with respect to the following 
Poisson brackets 



[Mi , Mj } ag = \ a k g x M k + x k ^ a m M m d m g 1 



3 



19 aid j 



i M * ' X 3 )ag = £ ^ 9 1 ^ 1 X k , { Xl , Xj } ag = . 

Reduction of the Poisson brackets {., .} ag to canonical Lie-Poisson brackets (|1.19j) and equivalence of 
the Veselova system to the Chaplygin ball without spinning are discussed in [7j. 
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5 Generalized Veselova system 

The Veselova system is an LR system on the Lie group G = SO(3) [T31I12] or on the Lie algebroid [S]. 
In " mathematical" equations of motion (|4.2[) were rewritten in the " physical" form 

Ilu = Ioj X LU + Xx , X = X X UJ , 

where I = A -1 being the inertia tensor of the rigid body and undetermined Lagrangian multiplier 

(I _1 :r, Iuj x uj) 



X 



(I-^z) 



was obtained from the mechanical reason stating an absence of the spinning around vector x. 

In this section we continue to study the equations of motion (|4.2p in their initial form [27J [29] . 
However, we will consider generalized Veselova system at 

(w,x)=b, del, (5.4) 

instead of (|4.ip following to [T2] ■ It means that the projection of the angular velocity vector a; on a; is 
a non-zero constant b. 

Let us introduce vector K instead of the angular momentum vector M 



K = M — 



(E - A)M, x 



(x,x) 

so that initial equations of motion (|4.2j) for the Veselova system read as 

K = K x lj , x — x x uj , where ui — AM , (5.5) 
whereas constraint (|5.4p takes the form of 

(x,K) =b. 

Below we will identify (x, K) with the Casimir element of the Poisson bracket. 

It is easy to see that the Jacobian of this transformation of coordinates is proportional to the 
Jacobi multiplier 



det 



d(x, M) 



g 2 (x) = a\x\ + cl2x\ + 03X3 . 



d(x,K) 

So, the transformation of coordinates M — > K is invertible where invariant measure [i ves (|4.3|) exists 
and where change of time exists, which allows us to integrate equations of motion in quadratures [271 . 

In variables x, K invariant measure of the equations of motion (15 ,5|) is expressed via old measure 
fj, ves and Jacobian 

fJiv = M-uesJ = 9 dxdK , g(x) = (x,Ax) . (5-6) 
It is easy to prove that integrals of motion 

Hi — (K, AgK) , H 2 = (K,K), H 3 = (x,x), H i = (x,K), (5.7) 

where 

A g = A - g~ 2 (E - A) (a; ® x) (E — A) , 
are in the involution with respect to the Poisson brackets 



3 



[Ki ,Kj} v = e ijk gK k + e ijk x k ^ K m d m g 



{ R i ) X J } v = £ ijk gXk, {Xi , Xj } v = . 



b 

9, 



(5.8) 



As usual (x, x) and (x, K) are the Casimir functions for these brackets and non-Hamiltonian equations 
of motion (|5.5p 

j t z = g- 1 (x){H,z} v , H=^, (5.9) 
become Hamiltonian equations after the change of time dt — > gdt. 



11 



Proposition 4 If [x, x) = 1 and parameters on of the Veselova system are related with parameters a-j 
of the Chaplygin system 

a,i = 1 — dai , (5.10) 

then the Poisson brackets {., .} v \5.8\) and invariant measure fi v \5.6\) coincide with the Poisson brackets 
{-,} g VS.8\) and invariant measure fi c h \3. 7[ ). 

In this case integrals of motion Hi |5. 7| ) for the Veselova system are expressed via integrals of 
motion Hi 13. 6\) for the Chaplygin ball 

H x =H 2 -dH x , H 2 = H 2 , H 3 = H 3 , H A = H A . (5.11) 

The proof is a straightforward. □ 

Consequently, the Veselova system the Poisson brackets {., .}„ are reduced to canonical Lie-Poisson 
brackets on e*(3) using the combination of (|3 . 11 1) and (|5.10[) 

Li = g- 1 
L 2 = g- 1 
L 3 = g- 1 

Inverse mapping is the combination of (|3.12l) and (I5.10|) as well. 

If we put c = 0, the combination of transformations (|5.9p and (|5.12l) relates the nonholonomic 
Veselova system with integrable Hamiltonian system on the sphere (|3.13p . The corresponding integrals 
of motion look like 

H 1 = H a - dHx , H 2 = H 2 . (5.13) 

So, nonholonomic Veselova system is a trajectory equivalent to the Hamiltonian system on the sphere 
(I3.13[) and, therefore, to the nonholonomic Chaplygin ball. These three dynamical systems have a 
common Lagrangian foliation and the common trajectories considered us unparametrized curves. 

Note that this equivalence explains the applicability of the Chaplygin transformation [10] to the 
Veselova system [12]. 

References 

[1] R. Abraham, J. Marsden: Foundations of Mechanics, (Second edition), Benjamin-Cummings, 
1978. 

[2] O. I. Bogoyavlenskii, Integrable cases of the dynamics of a rigid body, and integrable systems on 
the spheres Math. USSR-Izv., v.27, p. 203-218, 1986. 

[3] A.V. Bolsinov, B. Jovanovic, Complete involutive algebras of functions on cotangent bundles of 
homogeneous spaces. Math. Zeit., v. 246, p. 213-236, 2004. 

[4] A.V. Bolsinov, Complete commutative families of polynomials in Poisson-Lie algebras: A proof of 
the Mischenko-Fomenko conjecture, In book: Tensor and Vector Analysis, Vol. 26, Moscow State 
University, 2005, p. 87-109. 

[5] A.V. Bolsinov, B. Jovanovic, Magnetic flows on homogeneous spaces, Comment. Math. Hclv., 
v.83, p.679-700, 2008. 

[6] A.V. Borisov, I.S. Mamaev, The Chaplygin Problem of the Rolling Motion of a Ball Is Hamilto- 
nian, Math. Notes, 2001, vol. 70, no. 5, p. 720-723. 

[7] A.V. Borisov, I.S. Mamaev, Conservation Laws, Hierarchy of Dynamics and Explicit Integration 
of Nonholonomic Systems, Reg. Chaotic Dyn., v. 13(5), p. 443-490, 2008. 



K 2 - 



bx\ 
(x, x) 

bx 2 
(x,x) 

bx 3 
(x,x) 



V 

2 

l + f3 
V 



CX\ 



■ n -r x 2 
CX 2 



^2 



(5.12) 



• r 2 



v = h\x\ + a 2 x\ 



12 



A. V. Borisov, A. A. Kilin, I.S. Mamaev, Generalized Chaplygins transformation and explicit inte- 
gration of a system with a spherical support,. Rus. J. Nonlin. Dyn. , 2011, V.7, 2, p. 313-338. 

J. Cortes, M. de Leon, J.C. Marrero, E. Martinez, Nonholonomic Lagrangian systems on Lie 
algebroids, Discrete and Continuous Dynamical Systems, Series A, v. 24 (2), p. 213-271, 2009. 

S.A. Chaplygin, On a Ball's Rolling on a Horizontal Plane, Mathematical collection of the Moscow 
Mathematical Society, v. 24, p.139-168, 1903. 

B. Jovanovic, Hamiltonization and integrability of the Chaplygin sphere in R n , J. of Nonlinear 
Science, v.20, n.5, p.569-593, 2010. 

Yu. N. Fedorov, Two integrable nonholonomic systems in classical dynamics, Vestnik Moskov. 
Univ. Ser. I Mat. Mekh., 1989, no. 4, p. 38-41. 

Yu. N. Fedorov, B. Jovanovic, Hamiltonization of the generalized Veselova LR system, Regul. 
Chaotic Dyn., v. 14, p. 495-505, 2009. 

S. Hochgerner, Chaplygin systems associated to Cartan decompositions of semi-simple Lie algebras, 
Diff. Geom. App., v.28, n.4, p. 436-453, 2010. 

A. Lichnerowicz, Les varietes de Poisson et leurs algebres de Lie associeees, J. Diff. Geom., 1977, 
v.12, p.253-300. 

E. G. Kalnins, Separation of variables for Riemannian spaces of constant curvature, Addison- 
Wesley Educational Publishers Inc., U.S., 1986. 

C. -M. Marie, A property of conformally Hamiltonian vector fields; application to the Kepler 
problem, Preprint: arXiv:1011.573lV 2. 2010. 

V.G. Marikhin, V.V. Sokolov, On the reduction of the pair of Hamiltonians quadratic in momenta 
to canonic form and real partial separation of variables for the Clebsch top, Rus. J. Nonlin. Dyn., 
2008, vol. 4, no. 3, p. 313-322. 

T. Ohsawa, O.E. Fernandez, A.M. Bloch, D. V. Zenkov , Nonholonomic Hamilton- J acobi Theory 
via Chaplygin Hamiltonization J. Geometry and Physics, v. 61, p. 1263-1291, 2011. 

A. G. Reyman and M. A. Semenov-Tian-Shansky, Integrable systems: Group- Theoretical Ap- 
proach., Moskow-Izevsk, RCD, 2003. 

F. J. Turiel, Integrabilite dun tenseur de type (1,1) et structure symplectique du fibre cotangent, 
C.R. Acad. Sci., Paris, Ser. I, v. 301, p.923-925, 1985. 

A. V. Tsiganov, The Stdckel systems and algebraic curves, Journal of Mathematical Physics, v. 40, 
p.279-298, 1999. 

A. V. Tsiganov, Canonical transformations of the extended phase space, Toda lattices and Stdckel 
systems, Journal of Physics A, Math.Gen., v.33, p.4169-4182, 2000. 

A. V. Tsiganov, The Maupertuis principle and canonical transformations of the extended phase 
space, J. Nonlinear Math. Phys, v.8(l), p.157-182,2001. 

A. V. Tsiganov, Integrable Euler top and nonholonomic Chaplygin ball , Journal of Geometric 
Mechanics, v.3, n.3, p.337 - 362, 2011. 

A.V. Tsiganov, One invariant measure and different Poisson brackets for two non-holonomic 
systems, Regul. Chaotic Dyn., v.17(1), p. 72-96, 2012. 

L.E. Veselova, New cases of the integrability of the equations of motion of a rigid body in the 
presence of a nonholonomic constraint, in Geometry, Differential Equations and Mechanics, Gos. 
Univ., Mekh.-Mat. Fak., Moscow, p. 64-68, 1986. 

A. P. Veselov, Time change in integrable systems, Vestnik Moskov. Univ. Ser. I Mat. Mekh. 1987, 
no. 5, 25-29, 104. 



13 



[29] A. P. Veselov, L.E. Veselova, Integrable Nonholonomic Systems on Lie Groups, Mat. Notes, v. 44, 
no. 5, p. 810-819, 1988. 

[30] S. Wojciechowski, Integrable one-particle potentials related to the Neumann systems and the 
J acobi problem of geodesic motion on an ellipsoid, Phys. Lett. A, v. 107, p. 106-111, 1985. 



14 



